Evolutionary and developmental specialization of foveal cell types in the marmoset.

In primates, high-acuity vision is mediated by the fovea, a small specialized central region of the retina. The fovea, unique to the anthropoid lineage among mammals, undergoes notable neuronal morphological changes during postnatal maturation. However, the extent of cellular similarity across anthropoid foveas and the molecular underpinnings of foveal maturation remain unclear. Here, we used high-throughput single-cell RNA sequencing to profile retinal cells of the common marmoset (Callithrix jacchus), an early divergent in anthropoid evolution from humans, apes, and macaques. We generated atlases of the marmoset fovea and peripheral retina for both neonates and adults. Our comparative analysis revealed that marmosets share almost all their foveal types with both humans and macaques, highlighting a conserved cellular structure among primate foveas. Furthermore, by tracing the developmental trajectory of cell types in the foveal and peripheral retina, we found distinct maturation paths for each. In-depth analysis of gene expression differences demonstrated that cone photoreceptors and Müller glia (MG), among others, show the greatest molecular divergence between these two regions. Utilizing single-cell ATAC-seq and gene-regulatory network inference, we uncovered distinct transcriptional regulations differentiating foveal cones from their peripheral counterparts. Further analysis of predicted ligand-receptor interactions suggested a potential role for MG in supporting the maturation of foveal cones. Together, these results provide valuable insights into foveal development, structure, and evolution.

Ancient origin of the rod bipolar cell pathway in the vertebrate retina.

Vertebrates rely on rod photoreceptors for vision in low-light conditions. The specialized downstream circuit for rod signalling, called the primary rod pathway, is well characterized in mammals, but circuitry for rod signalling in non-mammals is largely unknown. Here we demonstrate that the mammalian primary rod pathway is conserved in zebrafish, which diverged from extant mammals ~400 million years ago. Using single-cell RNA sequencing, we identified two bipolar cell types in zebrafish that are related to mammalian rod bipolar cell (RBCs), the only bipolar type that directly carries rod signals from the outer to the inner retina in the primary rod pathway. By combining electrophysiology, histology and ultrastructural reconstruction of the zebrafish RBCs, we found that, similar to mammalian RBCs, both zebrafish RBC types connect with all rods in their dendritic territory and provide output largely onto amacrine cells. The wiring pattern of the amacrine cells postsynaptic to one RBC type is strikingly similar to that of mammalian RBCs and their amacrine partners, suggesting that the cell types and circuit design of the primary rod pathway emerged before the divergence of teleost fish and mammals. The second RBC type, which forms separate pathways, was either lost in mammals or emerged in fish.

Comprehensive single-cell atlas of the mouse retina.

Single-cell RNA sequencing (scRNA-seq) has advanced our understanding of cellular heterogeneity at the single-cell resolution by classifying and characterizing cell types in multiple tissues and species. While several mouse retinal scRNA-seq reference datasets have been published, each dataset either has a relatively small number of cells or is focused on specific cell classes, and thus is suboptimal for assessing gene expression patterns across all retina types at the same time. To establish a unified and comprehensive reference for the mouse retina, we first generated the largest retinal scRNA-seq dataset to date, comprising approximately 190,000 single cells from C57BL/6J mouse whole retinas. This dataset was generated through the targeted enrichment of rare population cells via antibody-based magnetic cell sorting. By integrating this new dataset with public datasets, we conducted an integrated analysis to construct the Mouse Retina Cell Atlas (MRCA) for wild-type mice, which encompasses over 330,000 single cells. The MRCA characterizes 12 major classes and 138 cell types. It captured consensus cell type characterization from public datasets and identified additional new cell types. To facilitate the public use of the MRCA, we have deposited it in CELLxGENE, UCSC Cell Browser, and the Broad Single Cell Portal for visualization and gene expression exploration. The comprehensive MRCA serves as an easy-to-use, one-stop data resource for the mouse retina communities.

Integrating genetic regulation and single-cell expression with GWAS prioritizes causal genes and cell types for glaucoma.

Primary open-angle glaucoma (POAG), characterized by retinal ganglion cell death, is a leading cause of irreversible blindness worldwide. However, its molecular and cellular causes are not well understood. Elevated intraocular pressure (IOP) is a major risk factor, but many patients have normal IOP. Colocalization and Mendelian randomization analysis of >240 POAG and IOP genome-wide association study (GWAS) loci and overlapping expression and splicing quantitative trait loci (e/sQTLs) in 49 GTEx tissues and retina prioritizes causal genes for 60% of loci. These genes are enriched in pathways implicated in extracellular matrix organization, cell adhesion, and vascular development. Analysis of single-nucleus RNA-seq of glaucoma-relevant eye tissues reveals that the POAG and IOP colocalizing genes and genome-wide associations are enriched in specific cell types in the aqueous outflow pathways, retina, optic nerve head, peripapillary sclera, and choroid. This study nominates IOP-dependent and independent regulatory mechanisms, genes, and cell types that may contribute to POAG pathogenesis.

Evolution of neuronal cell classes and types in the vertebrate retina.

The basic plan of the retina is conserved across vertebrates, yet species differ profoundly in their visual needs1. Retinal cell types may have evolved to accommodate these varied needs, but this has not been systematically studied. Here we generated and integrated single-cell transcriptomic atlases of the retina from 17 species: humans, two non-human primates, four rodents, three ungulates, opossum, ferret, tree shrew, a bird, a reptile, a teleost fish and a lamprey. We found high molecular conservation of the six retinal cell classes (photoreceptors, horizontal cells, bipolar cells, amacrine cells, retinal ganglion cells (RGCs) and Müller glia), with transcriptomic variation across species related to evolutionary distance. Major subclasses were also conserved, whereas variation among cell types within classes or subclasses was more pronounced. However, an integrative analysis revealed that numerous cell types are shared across species, based on conserved gene expression programmes that are likely to trace back to an early ancestral vertebrate. The degree of variation among cell types increased from the outer retina (photoreceptors) to the inner retina (RGCs), suggesting that evolution acts preferentially to shape the retinal output. Finally, we identified rodent orthologues of midget RGCs, which comprise more than 80% of RGCs in the human retina, subserve high-acuity vision, and were previously believed to be restricted to primates2. By contrast, the mouse orthologues have large receptive fields and comprise around 2% of mouse RGCs. Projections of both primate and mouse orthologous types are overrepresented in the thalamus, which supplies the primary visual cortex. We suggest that midget RGCs are not primate innovations, but are descendants of evolutionarily ancient types that decreased in size and increased in number as primates evolved, thereby facilitating high visual acuity and increased cortical processing of visual information.

Evolutionary and Developmental Specialization of Foveal Cell Types in the Marmoset.

In primates, high-acuity vision is mediated by the fovea, a small specialized central region of the retina. The fovea, unique to the anthropoid lineage among mammals, undergoes notable neuronal morphological changes during postnatal maturation. However, the extent of cellular similarity across anthropoid foveas and the molecular underpinnings of foveal maturation remain unclear. Here, we used high throughput single cell RNA sequencing to profile retinal cells of the common marmoset ( Callithrix jacchus ), an early divergent in anthropoid evolution from humans, apes, and macaques. We generated atlases of the marmoset fovea and peripheral retina for both neonates and adults. Our comparative analysis revealed that marmosets share almost all its foveal types with both humans and macaques, highlighting a conserved cellular structure among primate foveas. Furthermore, by tracing the developmental trajectory of cell types in the foveal and peripheral retina, we found distinct maturation paths for each. In-depth analysis of gene expression differences demonstrated that cone photoreceptors and Müller glia, among others, show the greatest molecular divergence between these two regions. Utilizing single-cell ATAC-seq and gene-regulatory network inference, we uncovered distinct transcriptional regulations differentiating foveal cones from their peripheral counterparts. Further analysis of predicted ligand-receptor interactions suggested a potential role for Müller glia in supporting the maturation of foveal cones. Together, these results provide valuable insights into foveal development, structure, and evolution. Significance statement: The sharpness of our eyesight hinges on a tiny retinal region known as the fovea. The fovea is pivotal for primate vision and is susceptible to diseases like age-related macular degeneration. We studied the fovea in the marmoset-a primate with ancient evolutionary ties. Our data illustrated the cellular and molecular composition of its fovea across different developmental ages. Our findings highlighted a profound cellular consistency among marmosets, humans, and macaques, emphasizing the value of marmosets in visual research and the study of visual diseases.

Integrated multi-omics single cell atlas of the human retina.

Single-cell sequencing has revolutionized the scale and resolution of molecular profiling of tissues and organs. Here, we present an integrated multimodal reference atlas of the most accessible portion of the mammalian central nervous system, the retina. We compiled around 2.4 million cells from 55 donors, including 1.4 million unpublished data points, to create a comprehensive human retina cell atlas (HRCA) of transcriptome and chromatin accessibility, unveiling over 110 types. Engaging the retina community, we annotated each cluster, refined the Cell Ontology for the retina, identified distinct marker genes, and characterized cis-regulatory elements and gene regulatory networks (GRNs) for these cell types. Our analysis uncovered intriguing differences in transcriptome, chromatin, and GRNs across cell types. In addition, we modeled changes in gene expression and chromatin openness across gender and age. This integrated atlas also enabled the fine-mapping of GWAS and eQTL variants. Accessible through interactive browsers, this multimodal cross-donor and cross-lab HRCA, can facilitate a better understanding of retinal function and pathology.

Suppressing DNMT3a Alleviates the Intrinsic Epigenetic Barrier for Optic Nerve Regeneration and Restores Vision in Adult Mice.

The limited regenerative potential of the optic nerve in adult mammals presents a major challenge for restoring vision after optic nerve trauma or disease. The mechanisms of this regenerative failure are not fully understood1,2. Here, through small-molecule and genetic screening for epigenetic modulators3, we identify DNA methyltransferase 3a (DNMT3a) as a potent inhibitor of axon regeneration in mouse and human retinal explants. Selective suppression of DNMT3a in retinal ganglion cells (RGCs) by gene targeting or delivery of shRNA leads to robust, full-length regeneration of RGC axons through the optic nerve and restoration of vision in adult mice after nerve crush injury. Genome-wide bisulfite and transcriptome profiling in combination with single nucleus RNA-sequencing of RGCs revealed selective DNA demethylation and reactivation of genetic programs supporting neuronal survival and axonal growth/regeneration by DNMT3a deficiency. This was accompanied by the suppression of gene networks associated with apoptosis and inflammation. Our results identify DNMT3a as the central orchestrator of an RGC-intrinsic mechanism that limits optic nerve regeneration. Suppressing DNMT3a expression in RGCs unlocks the epigenetic switch for optic nerve regeneration and presents a promising therapeutic avenue for effectively reversing vision loss resulted from optic nerve trauma or diseases.

Ancient origin of the rod bipolar cell pathway in the vertebrate retina.

Vertebrates rely on rod photoreceptors for vision in low-light conditions. Mammals have a specialized downstream circuit for rod signaling called the primary rod pathway, which comprises specific cell types and wiring patterns that are thought to be unique to this lineage. Thus, it has been long assumed that the primary rod pathway evolved in mammals. Here, we challenge this view by demonstrating that the mammalian primary rod pathway is conserved in zebrafish, which diverged from extant mammals ~400 million years ago. Using single-cell RNA-sequencing, we identified two bipolar cell (BC) types in zebrafish that are related to mammalian rod BCs (RBCs) of the primary rod pathway. By combining electrophysiology, histology, and ultrastructural reconstruction of the zebrafish RBCs, we found that, like mammalian RBCs, both zebrafish RBC types connect with all rods in their dendritic territory, and provide output largely onto amacrine cells. The wiring pattern of the amacrine cells post-synaptic to one RBC type is strikingly similar to that of mammalian RBCs, suggesting that the cell types and circuit design of the primary rod pathway have emerged before the divergence of teleost fish and amniotes. The second RBC type, which forms separate pathways, is either lost in mammals or emerged in fish.

Ancient origin of the rod bipolar cell pathway in the vertebrate retina.

Vertebrates rely on rod photoreceptors for vision in low-light conditions1. Mammals have a specialized downstream circuit for rod signaling called the primary rod pathway, which comprises specific cell types and wiring patterns that are thought to be unique to this lineage2-6. Thus, it has been long assumed that the primary rod pathway evolved in mammals3,5-7. Here, we challenge this view by demonstrating that the mammalian primary rod pathway is conserved in zebrafish, which diverged from extant mammals ~400 million years ago. Using single-cell RNA-sequencing, we identified two bipolar cell (BC) types in zebrafish that are related to mammalian rod BCs (RBCs) of the primary rod pathway. By combining electrophysiology, histology, and ultrastructural reconstruction of the zebrafish RBCs, we found that, like mammalian RBCs8, both zebrafish RBC types connect with all rods and red-cones in their dendritic territory, and provide output largely onto amacrine cells. The wiring pattern of the amacrine cells post-synaptic to one RBC type is strikingly similar to that of mammalian RBCs. This suggests that the cell types and circuit design of the primary rod pathway may have emerged before the divergence of teleost fish and amniotes (mammals, bird, reptiles). The second RBC type in zebrafish, which forms separate pathways from the first RBC type, is either lost in mammals or emerged in fish to serve yet unknown roles.

Transcriptomic analysis of the ocular posterior segment completes a cell atlas of the human eye.

Although the visual system extends through the brain, most vision loss originates from defects in the eye. Its central element is the neural retina, which senses light, processes visual signals, and transmits them to the rest of the brain through the optic nerve (ON). Surrounding the retina are numerous other structures, conventionally divided into anterior and posterior segments. Here, we used high-throughput single-nucleus RNA sequencing (snRNA-seq) to classify and characterize cells in six extraretinal components of the posterior segment: ON, optic nerve head (ONH), peripheral sclera, peripapillary sclera (PPS), choroid, and retinal pigment epithelium (RPE). Defects in each of these tissues are associated with blinding diseases-for example, glaucoma (ONH and PPS), optic neuritis (ON), retinitis pigmentosa (RPE), and age-related macular degeneration (RPE and choroid). From ~151,000 single nuclei, we identified 37 transcriptomically distinct cell types, including multiple types of astrocytes, oligodendrocytes, fibroblasts, and vascular endothelial cells. Our analyses revealed a differential distribution of many cell types among distinct structures. Together with our previous analyses of the anterior segment and retina, the data presented here complete a "Version 1" cell atlas of the human eye. We used this atlas to map the expression of >180 genes associated with the risk of developing glaucoma, which is known to involve ocular tissues in both anterior and posterior segments as well as the neural retina. Similar methods can be used to investigate numerous additional ocular diseases, many of which are currently untreatable.

Transcriptomic Analysis of the Ocular Posterior Segment Completes a Cell Atlas of the Human Eye.

Although the visual system extends through the brain, most vision loss originates from defects in the eye. Its central element is the neural retina, which senses light, processes visual signals, and transmits them to the rest of the brain through the optic nerve (ON). Surrounding the retina are numerous other structures, conventionally divided into anterior and posterior segments. Here we used high-throughput single nucleus RNA sequencing (snRNA-seq) to classify and characterize cells in the extraretinal components of the posterior segment: ON, optic nerve head (ONH), peripheral sclera, peripapillary sclera (PPS), choroid, and retinal pigment epithelium (RPE). Defects in each of these tissues are associated with blinding diseases - for example, glaucoma (ONH and PPS), optic neuritis (ON), retinitis pigmentosa (RPE), and age-related macular degeneration (RPE and choroid). From ∼151,000 single nuclei, we identified 37 transcriptomically distinct cell types, including multiple types of astrocytes, oligodendrocytes, fibroblasts, and vascular endothelial cells. Our analyses revealed a differential distribution of many cell types among distinct structures. Together with our previous analyses of the anterior segment and retina, the new data complete a "Version 1" cell atlas of the human eye. We used this atlas to map the expression of >180 genes associated with the risk of developing glaucoma, which is known to involve ocular tissues in both anterior and posterior segments as well as neural retina. Similar methods can be used to investigate numerous additional ocular diseases, many of which are currently untreatable.

Evolution of neuronal cell classes and types in the vertebrate retina.

The basic plan of the retina is conserved across vertebrates, yet species differ profoundly in their visual needs (Baden et al., 2020). One might expect that retinal cell types evolved to accommodate these varied needs, but this has not been systematically studied. Here, we generated and integrated single-cell transcriptomic atlases of the retina from 17 species: humans, two non-human primates, four rodents, three ungulates, opossum, ferret, tree shrew, a teleost fish, a bird, a reptile and a lamprey. Molecular conservation of the six retinal cell classes (photoreceptors, horizontal cells, bipolar cells, amacrine cells, retinal ganglion cells [RGCs] and Muller glia) is striking, with transcriptomic differences across species correlated with evolutionary distance. Major subclasses are also conserved, whereas variation among types within classes or subclasses is more pronounced. However, an integrative analysis revealed that numerous types are shared across species based on conserved gene expression programs that likely trace back to the common ancestor of jawed vertebrates. The degree of variation among types increases from the outer retina (photoreceptors) to the inner retina (RGCs), suggesting that evolution acts preferentially to shape the retinal output. Finally, we identified mammalian orthologs of midget RGCs, which comprise >80% of RGCs in the human retina, subserve high-acuity vision, and were believed to be primate-specific (Berson, 2008); in contrast, the mouse orthologs comprise <2% of mouse RGCs. Projections both primate and mouse orthologous types are overrepresented in the thalamus, which supplies the primary visual cortex. We suggest that midget RGCs are not primate innovations, but descendants of evolutionarily ancient types that decreased in size and increased in number as primates evolved, thereby facilitating high visual acuity and increased cortical processing of visual information.

Temporal single-cell atlas of non-neuronal retinal cells reveals dynamic, coordinated multicellular responses to central nervous system injury.

Non-neuronal cells are key to the complex cellular interplay that follows central nervous system insult. To understand this interplay, we generated a single-cell atlas of immune, glial and retinal pigment epithelial cells from adult mouse retina before and at multiple time points after axonal transection. We identified rare subsets in naive retina, including interferon (IFN)-response glia and border-associated macrophages, and delineated injury-induced changes in cell composition, expression programs and interactions. Computational analysis charted a three-phase multicellular inflammatory cascade after injury. In the early phase, retinal macroglia and microglia were reactivated, providing chemotactic signals concurrent with infiltration of CCR2+ monocytes from the circulation. These cells differentiated into macrophages in the intermediate phase, while an IFN-response program, likely driven by microglia-derived type I IFN, was activated across resident glia. The late phase indicated inflammatory resolution. Our findings provide a framework to decipher cellular circuitry, spatial relationships and molecular interactions following tissue injury.

Encoding of environmental illumination by primate melanopsin neurons.

Light regulates physiology, mood, and behavior through signals sent to the brain by intrinsically photosensitive retinal ganglion cells (ipRGCs). How primate ipRGCs sense light is unclear, as they are rare and challenging to target for electrophysiological recording. We developed a method of acute identification within the live, ex vivo retina. Using it, we found that ipRGCs of the macaque monkey are highly specialized to encode irradiance (the overall intensity of illumination) by blurring spatial, temporal, and chromatic features of the visual scene. We describe mechanisms at the molecular, cellular, and population scales that support irradiance encoding across orders-of-magnitude changes in light intensity. These mechanisms are conserved quantitatively across the ~70 million years of evolution that separate macaques from mice.

Vision-Dependent and -Independent Molecular Maturation of Mouse Retinal Ganglion Cells.

The development and connectivity of retinal ganglion cells (RGCs), the retina's sole output neurons, are patterned by activity-independent transcriptional programs and activity-dependent remodeling. To inventory the molecular correlates of these influences, we applied high-throughput single-cell RNA sequencing (scRNA-seq) to mouse RGCs at six embryonic and postnatal ages. We identified temporally regulated modules of genes that correlate with, and likely regulate, multiple phases of RGC development, ranging from differentiation and axon guidance to synaptic recognition and refinement. Some of these genes are expressed broadly while others, including key transcription factors and recognition molecules, are selectively expressed by one or a few of the 45 transcriptomically distinct types defined previously in adult mice. Next, we used these results as a foundation to analyze the transcriptomes of RGCs in mice lacking visual experience due to dark rearing from birth or to mutations that ablate either bipolar or photoreceptor cells. 98.5% of visually deprived (VD) RGCs could be unequivocally assigned to a single RGC type based on their transcriptional profiles, demonstrating that visual activity is dispensable for acquisition and maintenance of RGC type identity. However, visual deprivation significantly reduced the transcriptomic distinctions among RGC types, implying that activity is required for complete RGC maturation or maintenance. Consistent with this notion, transcriptomic alternations in VD RGCs significantly overlapped with gene modules found in developing RGCs. Our results provide a resource for mechanistic analyses of RGC differentiation and maturation, and for investigating the role of activity in these processes.

Cell atlas of the human ocular anterior segment: Tissue-specific and shared cell types.

The anterior segment of the eye consists of the cornea, iris, ciliary body, crystalline lens, and aqueous humor outflow pathways. Together, these tissues are essential for the proper functioning of the eye. Disorders of vision have been ascribed to defects in all of them; some disorders, including glaucoma and cataract, are among the most prevalent causes of blindness in the world. To characterize the cell types that compose these tissues, we generated an anterior segment cell atlas of the human eye using high-throughput single-nucleus RNA sequencing (snRNAseq). We profiled 195,248 nuclei from nondiseased anterior segment tissues of six human donors, identifying >60 cell types. Many of these cell types were discrete, whereas others, especially in the lens and cornea, formed continua corresponding to known developmental transitions that persist in adulthood. Having profiled each tissue separately, we performed an integrated analysis of the entire anterior segment, revealing that some cell types are unique to a single structure, whereas others are shared across tissues. The integrated cell atlas was then used to investigate cell type-specific expression patterns of more than 900 human ocular disease genes identified through either Mendelian inheritance patterns or genome-wide association studies.

Unified classification of mouse retinal ganglion cells using function, morphology, and gene expression.

Classification and characterization of neuronal types are critical for understanding their function and dysfunction. Neuronal classification schemes typically rely on measurements of electrophysiological, morphological, and molecular features, but aligning such datasets has been challenging. Here, we present a unified classification of mouse retinal ganglion cells (RGCs), the sole retinal output neurons. We use visually evoked responses to classify 1,859 mouse RGCs into 42 types. We also obtain morphological or transcriptomic data from subsets and use these measurements to align the functional classification to publicly available morphological and transcriptomic datasets. We create an online database that allows users to browse or download the data and to classify RGCs from their light responses using a machine learning algorithm. This work provides a resource for studies of RGCs, their upstream circuits in the retina, and their projections in the brain, and establishes a framework for future efforts in neuronal classification and open data distribution.

Overlapping transcriptional programs promote survival and axonal regeneration of injured retinal ganglion cells.

Injured neurons in the adult mammalian central nervous system often die and seldom regenerate axons. To uncover transcriptional pathways that could ameliorate these disappointing responses, we analyzed three interventions that increase survival and regeneration of mouse retinal ganglion cells (RGCs) following optic nerve crush (ONC) injury, albeit not to a clinically useful extent. We assessed gene expression in each of 46 RGC types by single-cell transcriptomics following ONC and treatment. We also compared RGCs that regenerated with those that survived but did not regenerate. Each intervention enhanced survival of most RGC types, but type-independent axon regeneration required manipulation of multiple pathways. Distinct computational methods converged on separate sets of genes selectively expressed by RGCs likely to be dying, surviving, or regenerating. Overexpression of genes associated with the regeneration program enhanced both survival and axon regeneration in vivo, indicating that mechanistic analysis can be used to identify novel therapeutic strategies.